Notes on SAW Tag Interrogation Techniques 


1. Introduction 

We consider the problem of interrogating a single SAW RFID tag with a known ID and 
known range in the presence of multiple interfering tags under the following 
assumptions: 

• The RF propagation environment is well approximated as a simple delay channel 
with geometric power-decay constant a > 2 . 

• The interfering tag IDs are unknown but well approximated as independent, 
identically distributed random samples from a probability distribution of tag ID 
waveforms with known second-order properties, and the tag of interest is drawn 
independently from the same distribution. 

• The ranges of the interfering tags are unknown but well approximated as 
independent, identically distributed realizations of a random variable p with a 
known probability distribution f p , and the tag ranges are independent of the tag 

ID waveforms. 

In particular, we model the tag waveforms as random impulse responses from a wide- 
sense-stationary, uncorrelated-scattering (WSSUS) fading channel with known 
bandwidth and scattering function. A brief discussion of the properties of such channels 
and the notation used to describe them in this document is given in the Appendix. 

Under these assumptions, we derive the expression for the output signal-to-noise ratio 
(SNR) for an arbitrary combination of transmitted interrogation signal and linear receiver 
filter. Based on this expression, we derive the optimal interrogator configuration (i.e., 
transmitted signal/receiver filter combination) in the two extreme noise/interference 
regimes, i.e., noise-limited and interference-limited, under the additional assumption that 
the coherence bandwidth of the tags is much smaller than the total tag bandwidth. Finally, 
we evaluate the performance of both optimal interrogators over a broad range of 
operating scenarios using both numerical simulation based on the assumed model and 
Monte Carlo simulation based on a small sample of measured tag waveforms. The 
performance evaluation results not only provide guidelines for proper interrogator design, 
but also provide some insight on the validity of the assumed signal model. 

It should be noted that the assumption that the impulse response of the tag of interest is 
known precisely implies that the temperature and range of the tag are also known 
precisely, which is generally not the case in practice. However, analyzing interrogator 
performance under this simplifying assumption is much more straightforward and still 
provides a great deal of insight into the nature of the problem. 

2. Basic Approach 

Let Hq (/) represent the known frequency response function (i.e., the Fourier transform 
of the impulse response) of the tag of interest, and let H k ( / ) for k = 1, 2, . . . , K, 



represent the unknown frequency response functions of the K interfering tags. We assume 
that all tag frequency response functions satisfy 

H k (/) = 0, V|/-/ 0 | > B/2, A: = 0, 1 , K , 

where B is the two-sided bandwidth of the tags and / 0 is an arbitrary center frequency. 
For the sake of notational simplicity, we will assume for the remainder of this document 
that we are operating with baseband-equivalent signals, so that f 0 = 0 . This involves no 
loss of generality as the results are all valid as stated for arbitrary center frequencies. 

If an arbitrary interrogating signal with frequency spectrum and energy £ is 

transmitted from the transmitting antenna, then the return signal from tag k = 0, 1, ..., K, 
received at the (co-located) receiving antenna is given by 

**(/)= 



where c is the signal propagation velocity and r k is the range of the tag from both the 
transmitting and receiving antennas measured in meters 1 . If the receiver implements a 
filter with frequency response H R (/) sampled at output time 1 = 0 , then the output from 

the receiver due only to the return signal X k (/) will be 


where £ T = J ^ | H T (/)| 2 df . Hence, the output from the receiver filter due to the signal 


of interest is given by 




H r(f) 


-i4nf r 0 /c 


and the aggregate interference from all of the interfering tags at the output of the receiver, 
is given by 


1 Note that the amplitude decay rate for the return signal is r a rather than r a ^ 2 since 

this is a two-way or cooperative radar channel. 
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Note that the aggregate interference represents a random variable with mean zero and 
variance 
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where cp(/-X) is the spaced-frequency correlation function (see Appendix) for the 
tags, and 

p(/-A,) = E jp -2a e -/2,r ^ - ^ p/c J 

represents a similar spaced-frequency correlation function for the multiplicative power 
decay and phase shift of the interfering tags. 

Finally, the output from the receiver filter due to a comp lex- valued additive white 
Gaussian noise (AWGN) process with power spectral density iV 0 and associated 

orthogonal increments frequency-domain process Z(/) is given by 

N = \% H B(f) Z W). 

which is a Gaussian random variable with mean zero and variance 


4=v 0 |_^ 2 |i/ s (/)| 2 rf/. 



Hence, if we let Y = S 0 +I + N represent the random output from the receiver filter in the 
presence of the signal of interest, additive interference, and AWGN, then it follows that 
the output SNR at the receiver is given by 
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( 1 ) 


Assuming that the aggregate interference is well approximated as Gaussian, the detection 
performance of the interrogator is determined by this quantity. Such an assumption is 
probably only justified for a large number of interfering tags, but maximizing the output 
SNR remains a worthwhile optimization criterion in any case. 

As a rule, Expression (1) can be simplified quite a bit by making the following additional 
assumptions, which should be valid approximately for tags with good anti-collision 
properties and appropriate choices of H T (/) and H R {/)'■ 


1. (p(o) = 0V|u| > A/2 , where A « B . Recall that A represents the coherence 
bandwidth for the set of tag waveforms, 


2. H R (f)»H R (X),H T (f)~H T (X)V\f-X\<A, 

3. H k (f)*H k (X)\/\f-X\<A,k = 0,l,...,K. 


Under these assumptions, the expression for the variance becomes 
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which is a considerable simplification, and the expression for the output SNR becomes 
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where iVy = J ^(p(u)p(u)£/u. Note that as K£Nj/N 0 — »oo, the interference term in 
the denominator dominates, and this expression simplifies to 
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On the other hand, as K£Nj /N 0 — > 0 , the AWGN term dominates, and the expression 
simplifies to 
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3. Optimal Interrogation 

Note that in the general case, the choice of the transmitted signal H T (/) and the 

receiver filter H R [f ) that maximize the output SNR given by Expression (1) must be 

determined numerically by solving a fairly complex variational problem. Hence, in 
general, the optimal interrogator design and the corresponding detection performance is 
rather difficult to determine. However, for practical purposes, the problem can often be 
reduced to either the interference-limited case or the noise-limited case 2 . In each of the 
two limiting cases, the optimal interrogator design and corresponding optimal output 
SNR are easy to determine. We consider each case separately below. 


2 The problem reduces to the interference-limited case as the product of the number of 
interfering tags and the power of the transmitted signal goes to infinity. As the same 
product goes to zero, the problem reduces to the noise-limited case. 



3.1. Noise-Limited Case 

For this case, the output SNR is given by Expression (3) and satisfies 


SNR N =r„- 
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It follows from the Cauchy-Schwarz (C-S) inequality that equality is achieved in this 
expression if and only if 

Hr (/) = (3^2# o (f)e i4n/r ° /c , 

for an arbitrary complex constant P * 0 . Furthermore, in this case, we have 

sNR »= r " a ^C|^r i// »( / )| dfir °- 2a ^O H °( f t df - 

where equality is achieved if and only if 

n T {f)=yH o(/). 

for an arbitrary constant y ^ 0 . Hence, without loss of generality, the optimal 
interrogating signal satisfies 




the optimal receiver filter satisfies 


H R (f)= 

where £ 0 = [ff df and p ^ 0 is arbitrary, and the optimal output SNR is given 
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Finally, averaging over the constellation of random tag waveforms gives 


sm N = , 0 - 2 “ {J_; /2 |/f 0 (/)| df } = r 0 - 2a %(o)— . 



Notice that in this case, the optimal transmitted signal is just the time-reverse of the 
impulse response of the interrogating signal for the tag of interest. This is not surprising 
given the “well-known” optimality properties of time-reversal processing [1,2]. 

3.2. Interference-Limited Case 


For this case, the output SNR is given by Expression (2) and satisfies 
\\ B 'L H R (f) H T (f)H 0 (/K' 4 ' / ' i)/c # 2 
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Again, it follows from the C-S inequality that equality is achieved in this expression if 
and only if 

H R {f)H T (f) = W 0 _ 


for an arbitrary complex constant p ^ 0 . Hence, one possible choice for the optimal 
interrogating signal satisfies 


»r(/) = 


j I. l/b-B/2, 
jo, |/|>£/2, 


the corresponding optimal receiver filter satisfies 

H R (f) = pH 0 (f)e i4nfr °' c , 

where P ^ 0 is arbitrary, and the optimal output SNR is given by 

SNR '=^Q H °( f f d/ - 

Finally, averaging over the constellation of random tag waveforms gives 

Notice that in this case, the optimal transmitted signal is just a band-limited impulse 
function, which can be approximated using either a stepped CW signal or an FM chirp. 

4. Performance Evaluation 

To compare the performance of the two interrogator configurations derived in Sections 
3.1 and 3.2 above (which would be optimal in the noise-limited case and the interference- 
limited case, respectively, if our assumptions were satisfied), we computed the value of 
the output SNR for each configuration both numerically and using Monte Carlo 
simulation over a wide range of operational scenarios. All perfonnance results were 
derived from a collection of 23 measured tag frequency response vectors that were 
generated by sampling the actual frequency response of 23 SAW RFID tags chosen 
randomly from a homogenous population of tags provided by a single manufacturer. 



Numerical performance results were derived by evaluating Expression (1) numerically 
using a sampled version of cp(/-X) estimated directly from the population of tag 

frequency response vectors. Simulated performance was derived by generating random 
realizations of the random variable 7 = S 0 + 1 + N that underlies Expression ( 1 ), directly 

estimating the quantities E{7] and Var ]7j , and computing an estimate of output SNR 

as E{7}/Var{7}. 

For the performance results presented in the following two sections, all tags (i.e., the tag 
of interest as well as all interfering tags) were assumed to be at the same range, 
represented by the variable Ro, propagation was assumed to be in free space (i.e., a = 2), 
the tag processing bandwidth was fixed at 80 MHz (indicated in the figure titles by the 
value of BW Factor =2), and the tag antenna temperature was assumed to be 300° Kelvin. 

The assumed gain of the transmitting antenna is represented by the variable G T , the gain 
of the receiving antenna is represented by the variable Gr, and the output SNR curves in 
the figures correspond to 21 different transmitted power levels uniformly spaced at 1 
dBm increments from 10 dBm to 30 dBm with a fixed receiver processing bandwidth 
(i.e., IF bandwidth) of 15 kHz. Finally, all tag frequency response waveforms were scaled 
to correspond to an insertion loss of approximately 30 dB in total tag input vs. output 
energy. It follows that for all figures, the output SNR value for the interference-optimal 
(i.e., band-limited impulse) interrogator with no interfering tags corresponds exactly to 
the value of SNR at the input to the receiver (i.e., E b /N 0 ) computed using a traditional 
link-budget calculation. 

4.1. Numerical Performance Results. 

For all figures in this section, the 21 curves corresponding to the noise-optimal (i.e., time- 
reversal) interrogator are approximately 4 dB above the corresponding interference- 
optimal interrogator curves when no interfering tags are present, but drop below the 
interference-optimal curves in many cases with five interfering tags. If the curves were 
extended to 100 interfering tags (excluded due to space limitations), both the noise- 
optimal and the interference-optimal curves would be essentially independent of 
transmitted power near the 100-tag extreme (i.e., only one curve for each instead of 21) 
and the noise-only curves would all fall below the corresponding interference-optimal 
curves at some point in the 0-100 interfering-tag range. This confirms that the noise- 
optimal interrogator consistently outperforms the interference-optimal interrogator in the 
noise-limited regime and vice-versa in the interference-limited regime. Interestingly, in 
many cases, the interference-limited regime is reached with only one interfering tag 
present. 
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Numerical Results 
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4.2. Numerical vs. Simulated Performance Results. 


In this section, we compare the numerical performance results with the simulated 
performance results as an indication of how well the assumed statistical tag-response 
model captures the behavior of the small number of measured tag responses in our 
sample. In all cases, the modeled behavior was nearly identical to the simulated behavior 
in the absence of interfering tags (as one would definitely expect), but the modeled output 
SNR always slightly exceeded the simulated SNR when interfering tags were present. 
Nevertheless, considering the very small size of our measured tag sample, the results 
were in remarkably good agreement throughout the range of 0-5 interfering tags. 


Numerical vs Simulated Results 
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Numerical vs Simulated Results 
Interference-Optimal Interrogators 
R0=30 m, Alpha=2, GT=3 dBi, GR=3 dBi, BWFactoi=2 
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Appendix. Background on Wide-Sense-Stationary Uncorrelated-Scattering Fading 
Channels 

In this appendix, we review the basic concepts and notation for the wide-sense-stationary, 
uncorrelated-scattering (WSSUS) statistical model of multipath fading channels that is 
utilized in the previous sections of this document. For a much more thorough discussion, 
see, for example, [3]. 

We assume that the channel under consideration has a time-varying impulse response 
given by h(x\t) , where t represents the time at which the output is measured and t 

represents the delay of the measurement. That is, h(x;t) represents the output of the 

channel at time t due to an impulse transmitted at time t-x. It is further assumed that 
h(x;t) is a zero-mean, complex random process that is causal, stationary in t, and 
satisfies the uncorrelated scattering assumption; i.e., 

E{/r(x;f)} =0 Vf,T, 

and 


E{Mx,;, 1 )r(x 2; , 2 )}=H T,; ''" 2)8(Tl ' T2) ’ T, -°' 


0, otherwise. 


where an overbar indicates complex conjugation, and 8(x) represents the Dirac delta 

function. Note that if we want to avoid 8 functions, we can express everything in the 
delay frequency domain. That is, if we define the Fourier transform of h(x\t) with 

respect to x as 

H(f-,t) = f_^h(x-t)e~ i2nfx dx, 


then we get 


E {jy ( Xi ; ^ ) ^ ( A ;^ 2 )} = E { /i ( x, ;^ ) A ( x 2 ; / 2 )} £ /x 1 z/t 2 

= j7<p( T f h -t 1 )e~ ,2m ' i<A ~ fl) dx i 
= <p(/i ~ fi>h -h)> 


where (p represents the so-called spaced-time, spaced-frequency correlation function. 
Note that cp is only a function of the difference in both time and frequency; that is, the 
random process H [f;t ) is stationary in both variables. If the explicit output time is not 
of interest, the t is frequently dropped in favor of the short-hand notation H (f ). 
Similarly, when q -t 2 =0, the notation <p(/i — / 2 ;0) is usually replaced by <p(/i ~/ 2 ), 
which is then referred to simply as the spaced frequency correlation function. 



Note also that cp is often described in terms of its 2 -D Fourier transform, which is given 
by 


5 (x;X) = O" e-‘ 2 ,ai dfd, 

The function is generally referred to as the scattering function of the channel and 

represents the power spectral density of the two-dimensional stationary random process 
//(/;/). Finally, it should be noted that we make the assumption of zero mean primarily 

for simplicity in exposition, but it is realistic in many cases, particularly if the channel 
has no line-of-sight component. In any case, extending the results to account for nonzero 
mean is straightforward. 

Now, if we transmit a signal .?(/) over such a channel, then, ignoring additive noise, the 
output is a zero-mean, complex random process of the form 

K0=O( T ;'M'- T y x ’ 

The autocorrelation function of this process is given by 

E{r(h)F(r 2 )} = Ej \j^(x 2 -,t 2 )j(t 2 -T 2 )dx 2 

= j:^i ; h-, 2 ).(h-xiWr 2 -xiK, 

Note that this autocorrelation function implies that the output is generally a nonstationary 
process unless the signal has some special properties. For example, if s(t\ is itself a zero 
mean, stationary random process that is independent of the channel process, then 

E {'• (h )~{h )} = 9 ( Ti ; t x - 1 2 ) E {5 (ft - Ti ) J( t 2 - Tj )} dx x 

= <P,('i ~h\ '2)^1 > 

where (h - ) represents the autocorrelation function of the signal. 
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